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Abstract 

If gaugino condensations occur in the early universe, domain walls are produced 
as a result of the spontaneous breaking of a discrete R symmetry. Those domain 
walls eventually annihilate with one another, producing the gravitational waves. We 
show that the gravitational waves can be a probe for measuring the gravitino mass, 
if the constant term in the superpotential is the relevant source of the discrete R 
symmetry breaking. 



1 Introduction 



The presence of numerous gauge symmetries is quite a mediocre phenomenon in the string 
landscape pQ and it is even unlikely that our universe possesses only the standard-model 
(SM) gauge interactions. Thus, we postulate, through this paper, that we have many 
(hidden) gauge bosons besides those in the SM. 

Each gauge boson is always accompanied by its fermionic partner called a gaugino in 
supersymmetry (SUSY). Some of the gauginos may condense in the vacuum by the strong 
forces of the corresponding gauge interactions. Indeed, the gaugino condensation is known 
to occur in e.g. a SU(N) pure SUSY Yang-Mills theory. In the theory, a continuous R 
symmetry is explicitly broken to a discrete Z 2N subgroup by instantons [2j. The gaugino 
condensation spontaneously breaks Z 2 n further down to Z 2 . Accordingly, domain walls 
are formed at epochs of the condensations in the early universe (3j HJ [5], EJ [7j. If the 
gaugino condensation occurs before the inflation, the domain walls are diluted away by 
the subsequent exponential expansion, and there is no cosmological remnant of the walls in 
the observable universe. However, if the gaugino condensation occurs after the inflation, 
the domain walls may play an important role in cosmology, since their energy density 
decreases much more slowly than ordinary matters. In fact, if the domain walls are 
stable, they will dominate the energy density of the universe in the end, which makes 
the universe extremely inhomogeneous. To avoid such cosmological disaster, the domain 
walls must eventually disappear. Fortunately, we know that the discrete Z 2 n symmetry 
is not an exact symmetry and hence the domain walls are not stable. 

In this paper, we discuss the dynamics of the domain walls formed by the gaugino 
condensation, assuming that the SU(N) gauge sector does not have sizable couplings to 
the other sectors. We will briefly discuss the case that such couplings to the SUSY break- 
ing sector are present later. Then the source of the breaking of the R symmetry relevant 
for the SU(N) sector comes only from the constant term wq in the superpotential W. 
The constant term preserves only the Z 2 symmetry, and it is needed to cancel the cosmo- 
logical constant in supergravity (SUGRA). We find that the walls collapse at H ~ 7773/2, 
where H is the expansion rate of the universe and 777,3/2 the gravitino mass. Interestingly, 
gravitational waves [S] are likely produced during the violent annihilation processes of 
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the domain walls |9J. The typical frequency of the gravitational waves we observe today 
depends on the gravitino mass, while the intensity is sensitive to the gaugino condensa- 
tion scale, A. Moreover, thermal history after inflation affects both the frequency and 
the intensity. We point out in this paper that such stochastic gravitational waves may be 
detectable by the forthcoming experiments such as Advanced LIGO [TU1 HI). LCGT [12], 
LISA [13], BBO [H] and DECIGO [15], for wide ranges of the gravitino mass and the 
reheating temperature. 

2 Gaugino condensation and domain walls 

We briefly explain the gaugino condensation and the subsequently formed domain walls. 
Let us consider a pure SU(N) SUSY Yang-Mills theory 0. Then the gauge interactions 
become strong at a dynamical scale A, and the gaugino condensation occurs: 

(X a X a ) k = -32n 2 e 2nik/N A 3 , (1) 

where k = 1,...,N, and there are N degenerate vacua. The effective superpotential 
below the dynamical scale is given by 

W GC = NA 3 e 2nik/N . (2) 

One can derive the forms of the above expressions using the symmetry and the holomorphy 
of the superpotential. 

We assume that the universe undergoes an exponential expansion called inflation. 
After inflation ends, the inflaton continues to dominate the energy density of the universe 
until it reheats the universe by the decay. Even before the reheating, there is a background 
thermal plasma. The maximal temperature of the background plasma can be expressed 
in terms of the reheating temperature Tr and the Hubble parameter during inflation H in { 
as [El EE] 

T^~(T 2 R H- mi M P f , (3) 

1 One may end up with a pure SUSY Yang-Mills theory at a low-energy scale, after integrating out 
heavy flavors. 
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where Mp = 2.4 x 10 18 GeV is the reduced Planck mass. We require that the maximal 
temperature of the background plasma is higher than the dynamical scale: 

T m > A. (4) 

Then, the non-zero gaugino condensate is formed after inflation, when the cosmic tem- 
perature becomes comparable to A 0. When gaugino condensation occurs, one of the N 
vacua is selected in each domain with the size of 0(£), where £ denotes the correlation 
length, £ ~ T~ x . Since the correlation length is much smaller than the horizon scale 
at the formation, there appear domain walls that connect the distinct vacua populated 
in each domain. The properties of the domain walls have been extensively studied (see 
Refs. [6], HE])- In particular, the tension of the wall connecting the adjacent vacua is 
analytically given by [6] 



a = 2NA 3 



e 2m/N 



(5) 



In the limit of large TV, a approaches 47rA 3 . In the following analysis we adopt a = 47rA 3 
as a reference value. 

How does the domain-wall network evolve? The numerical [191 Ell [21] and analyt- 
ical [22] results obtained so far show that the domain walls reach a scaling regime in 
simple models. We therefore assume that the evolution of the domain walls is given by 
the scaling law B That is, an averaged number of the walls per horizon remains the same 
in the evolution of the universe. This should be the case if the annihilation processes of 
the domain walls are efficient. The energy density of the domain walls, Pdwi obeying the 
scaling law is given by 

Pdw — naH, (6) 

where H is the Hubble parameter, n denotes an averaged number of the walls per horizon, 
and we expect n ~ 0(1). Note that the energy density of the domain walls decreases much 



2 Here and in what follows we assume that the SU (N) gauge sector is in thermal equilibrium with a 
temperature comparable to that of the SM sector. If the SU(N) gauginos are not in thermal equilibrium, 
gaugino condensation occurs when the Hubble parameter, not the temperature, becomes comparable to 
A. For the domain walls to be formed, H m f ^ A must be satisfied. Even in this case, if the domain-wall 
network reaches the scaling regime, our arguments remain essentially unchanged. 

3 We assume that there is no light particles in plasma which strongly couple to the domain walls. 
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more slowly than the radiation and ordinary matters. Therefore, if the domain walls are 
stable, they come to dominate the energy density of the universe when 

H ~ 4h- (7) 

Once the domain walls dominate the energy density of the universe, the universe will 
become intolerably inhomogeneous and anisotropic [5J. Therefore the domain walls must 
disappear before it starts to dominate the universe. 

One of the easiest ways to induce the decay of the domain walls is to introduce a 
bias and lift the degeneracy of the vacua. That is, if the spontaneously broken discrete 
symmetry is only approximate, the false vacua will decay into the true vacuum in the 
end, and the domain walls will disappear. Such biased domain walls was first studied 
by Vilenkin [5]. Here we summarize the result. Let e denote the bias between the false 
and true vacua, (e has mass dimension 4.) Then the domain walls annihilate with one 
another and disappear when the tension and the pressure on the wall due to the bias 
become comparable. The Hubble parameter at the decay is 

-f^decay ~ ~ • (8) 

Since the decay must be complete before the walls come to dominate the universe, there 
is a lower bound on the bias: 

On the other hand, if the bias is too strong, the domain walls are not formed, or even if 
formed, the walls might disappear before they reach the scaling law. We require that the 
bias is much smaller than ppw when they reach the scaling law soon after the formation. 
That is, 

oH f > e (10) 



with 



A 2 

for T R > A 



M K S , (11) 

where Hf denotes the Hubble parameter at the domain-wall formation. 
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Our concern here is if such a bias exists in the case of the gaugino condensation. 
Fortunately, there is indeed a source for the explicit breaking of the discrete Z 2 n symmetry. 
In order to cancel the cosmological constant, we need to introduce a constant Wq in the 
superpotential. The constant explicitly breaks Z 2 n down to Z 2 , and so, the degeneracy 
of the vacua is lifted. To see this, we write down the scalar potential in SUGRA using 
the effective superpotential W = Wqc + w o B, 

V ~ -3jV^A 3 e Jhri */ JV + h.c.. (12) 
Mp v ' 

Thus, depending on the relative phase between the gaugino condensate and w , the scalar 
potentials take different values in the distinct N vacua. The typical size of the bias 
between the adjacent vacua is given by 

e ~ am 3/2 A 3 , (13) 

where we have used \wq\ = m 3 / 2 Mp, and a = O(10) is a numerical coefficient. Throughout 
this paper, we assume that wo is the dominant source of the bias. Substituting (jSJ) and 
(TT3"1) into dSJ), one can see that the Hubble parameter at the decay is given by 

#decay ~ ™3/2- (14) 

Thus, as long as the constant wq is the main source of the bias, the timing of the domain- 
wall decay is determined by the gravitino mass. The constraints on the bias, (jHJ) and (IIP]) , 
read 

A^ in « A 3 < 0.1m 3/2 M P (15) 

with 



m 3/2 M P for Tr> A 

Amin e v j . 16 

{m 3/2 MpT 2 ) 4 for T R < A 



As we see in the next section, the intensity of the gravitational waves increases as A 
becomes large. When A takes the largest allowed value, the energy of the domain walls 
becomes comparable to the total energy of the universe at the decay. 



4 The non-zero value of Wq produces the gaugino mass, m\ = (3A^g 2 /167r 2 )m 3 / 2 , by the anomaly 
mediation [23] . where g is the gauge coupling of SU(N). One can also derive the effective potential, (|12p. 
from the gaugino mass term with the gaugino condensation [TB] , 
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The energy density stored in the domain walls must go into light degrees of freedom. 
In this paper we simply assume that the energy quickly goes into another hidden massless 
U(l) gauge bosons through loop diagrams of some heavy vector-like quarks charged under 
both gauge symmetries □. Since we assume that the domain walls decay before they come 
to dominate the universe, the energy of the hidden radiation is subdominant compared 
to the visible one. In Sec. HI we will discuss a possibility to observe the effect of the 
hidden radiation. Even if the domain walls decay into the visible particles, the following 
arguments on the gravitational waves remain intact, as long as the decay proceeds fast 
enough. 

Lastly let us summarize our picture. The domain walls are formed when the gauginos 
condense, i.e., T ~ A. Those walls evolve according to the scaling law ([6]) until H ~ 
When H ~ m 3 ^, the domain walls experience the pressure due to the bias, and they are 
pulled toward one another. After violent collisions of the domain walls, the false vacua 
decay into the true vacuum, and the domain walls disappear. The energy stored inside 
the domain walls and the false vacua will be released into light degrees of freedom. Most 
importantly, since the walls are spatially extended two dimensional objects, the collisions 
likely occur in such a way widely deviated from the spherical symmetry, leading to the 
production of the gravitational waves. Thus the domain walls of the gaugino condensate 
may have left their traces in the gravitational waves, which will be discussed in the next 
section. 

3 Gravitational waves 

We discuss here the properties of the gravitational waves produced from the decay of the 
domain walls. To detect the gravitational waves by the forthcoming experiments, the 
frequency and the intensity must be in certain ranges. 

Let us first consider the frequency of the gravitational waves. The gravitational waves 
are expected to peak at the frequency corresponding to a typical physical scale when they 
are produced. In our case, this frequency corresponds to the horizon scale at the decay, 



5 We assume that the effective couplings between the U(l) and SU(N) sectors, obtained after the 
integration of the heavy quarks, are so weak that the domain- wall dynamics is not affected. 
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i.e., /* ~ m 3 / 2 . The frequency is red-shifted due to the subsequent cosmic expansion, and 
the frequency we observe today, f , is much smaller than 

/o = — /*, (17) 
a 

where ao and a* are the scale factors at present and at the formation of the gravitational 
waves, respectively. The amount of the red-shift, or equivalently, a^/a , depends on the 
thermal history of the universe. 

If the reheating of the universe is completed before H ~ m 3 / 2 , the universe is radiation- 
dominated afterwards until the matter-radiation equality. The ratio of the scale factors 
is estimated by 

^ ~ 6 x 1(T 14 f-^V 5 ( , (18) 
a V200/ VIGeV/ ' V ; 

where T* denotes the cosmic temperature, and g* counts the relativistic degrees of freedom, 
both evaluated at a = a*. Thus the frequency at present, /o, is directly related to the 
gravitino mass as 

/o - 0.1Hz f^VV-^*. (19) 
J V200/ VlkeV/ V ; 

On the other hand, if the reheating is not completed when H ~ m 3 / 2 , the universe is 

matter-dominated until the reheating. Then the amount of the red-shift is modified as 

The present frequency becomes 

(iSvH^)'- «-) 

Surprisingly, the peak frequencies of the gravitational waves coincide with those covered 
by various ongoing and planned experiments, for wide ranges of the gravitino mass and 
the reheating temperature. We will come to this point later. 

The intensity of the gravitational waves decreases as the universe expands, since the 
amplitude is red-shifted for a mode with the wavelength less than the horizon size. In 
order to characterize the intensity, one may use a dimensionless quantity f2 gw (/) defined 

^ - yM- (22) 
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where p gw is the energy density of the gravitational waves, p c the critical energy density, 
and / the frequency. 

Let us estimate the magnitude of f2 gw (/) (see [8]). The energy of the gravitational 
waves in a horizon at H = m 3 / 2 is estimated by 

M 2 DW 2ttA 6 
E gw ~ G — — ~ —3 — —j, (23) 

where G = l/(8irMp) is the Newton constant, Mow the energy stored in the domain 
walls, and R* the typical spatial scale of the energy distribution. In the second equality 
in (|23|) . we have used Mp>w ~ cr/my 2 and R* ~ l/m 3 / 2 . If the reheating is completed 
before H = m 3 / 2 , we have 

^(fo)h 2 ~ 10- 7 f^V"f T )\ (24) 

where fo is given by (Fl9|) . and /i is the present Hubble parameter in units of 100 km/ s/Mpc. 
For A close to its upper bound given by (fT5]) . therefore, the energy density of the grav- 
itational wave amounts to one percent of the cosmic microwave background photons! If 
the reheating is not completed at H = m 3 / 2 , the gravitational waves are diluted by the 
entropy production of the inflaton. The density parameter, f2gw, is then given by 

T R y(m 3/2 yl( A 3 \ 2 
lOOGeV/ VlkeV/ V0.1m 3/2 Mj,/ ' v ; 

where fo is given by ( |2TI) . Note that the intensity of the gravitational waves is very 
sensitive to the dynamical scale, A. In order to detect the gravitational waves from the 
domain-wall decay, A must be close to its upper bound. 

We have made a rough estimate of the intensity of the gravitational waves in the 
above. Further detailed analyses on the domain-wall dynamics may reveal the presence of 
additional suppression factors. However, there are at least two reasons to believe that the 
gravitational waves produced by the domain-wall dynamics are not extremely suppressed. 
One is that the collisions among the domain walls are strongly first order phase transition, 
since they are decoupled from the thermal bath. The other is that the domain-walls are 
spatially extended two dimensional objects, which exhibits a peculiar dependence on the 
spatial scale. To see this more explicitly, let us consider the dependence of f2 gw on the 
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spatial scale, R. Using M^w oc R 2 , we have E gw oc R 3 , leading to f2 gw oc R°\ This is 
very different from the ordinary non-relativistic matters, which would lead to Q gw oc R 2 . 
Therefore, the resultant intensity of the gravitational waves is independent of the spatial 
scale, which makes the above rough estimate more robust 0. Thus, the domain-wall 
dynamics may be an ideal source for generating the gravitational waves in the early 
universe. 

Here let us briefly mention the sensitivities of the ongoing and planned experiments on 
the gravitational waves. One of the ground-based experiments, LIGO [24], is in operation 
and it is sensitive to the frequency between O(10) Hz and 10 4 Hz. The latest upper 
bound is fl^h 2 < 6.5 x 10~ 5 around 100 Hz [25], and an upgrade of the experiment, 
Advanced LIGO, would reach sensitivities of O(10~ 9 ). The sensitivity of LCGT would 
be more or less similar to that of Advanced LIGO. There are also planned space-borne 
interferometers such as LISA, BBO and DECIGO. LISA is sensitive to the band of (0.03 — 
0.1) mHz^S / ~0.1Hz, and it can reach Q gw h 2 < 10~ 12 at f = 1 mHz. Moreover, BBO 
and DECIGO will cover 10 mHz ^ /o ~ 10 2 Hz with much better sensitivity. 

In Fig.[TJ we show the contours of the frequencies: /o = 10~ 4 , 10~ 2 , 1, 10 2 , and 10 4 Hz, 
and the (maximal) intensities: Q gw h 2 = 10~ 22 , 10~ 17 , 10~ 12 and 10~ 7 , in the parameter 
space of the gravitino mass and the reheating temperature. Here we have set A to be the 
largest value allowed by (ITBT) . In the shaded triangle region, Q gw h 2 is saturated to 10~ 7 , 
since the reheating is complete before the domain- wall decay (see ( fl9l) and ( 1241) ). Also 
we schematically show the upper bounds on Tr from the gravitino problem [261 123 [28J 
EE EQl EH [321 [331 EU (solid (blue) line). The dotted (purple) line shows the constraint 
obtained by (gj) and H mi < 10~ 4 M P [35], and the region below the line is excluded B 

One can see that, for wide ranges of 7713/2 and Tr, the peak frequencies of the grav- 
itational waves correspond to those covered by the ongoing and planned experiments. 
Compared with the sensitivities of the experiments mentioned above, the gravitational 
waves from the domain- wall decay may be detectable, even by Advanced LIGO. Once it 
is detected, we will be able to obtain information on both the gravitino mass and the 
reheating temperature. In particular, for m 3 / 2 ^ O(10) eV, the measured frequency can 

6 The peak frequency scales as though. 

7 If the SU(N) sector is not thermalized, the constraint of the dotted (purple) line is not applied. 
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m 3 /2[GeV] 

Figure 1: Contours of the frequency (thick-dotted, red) and the maximal intensity 
(thin-dashed, green) of the gravitational waves from the domain-wall decay. A is set 
to be the largest value allowed by ffT5l) . The solid (blue) line schematically shows upper 
bounds on the reheating temperature from the gravitino problem. The dotted (pur- 
ple) line shows the bound coming from gD and H ini < 10" 4 M P , and the region below 
the dotted (purple) line is excluded. Note that the bound is not applied if the SU (N) 
gauginos are not in thermal equilibrium. LIGO and LISA are sensitive to the bands of 
O(10) Hz < f < 10 4 Hz and 0.1 mHz < /„ < 0.1 Hz, respectively, while BBO and DECIGO 
will cover 10 mHz ^5 /o ~ 10 2 Hz with much better sensitivity. 

be directly related to the gravitino mass by ffT9]) . 

4 Discussion and conclusions 

We have so far focused on the gravitational waves produced by the domain-wall dynamics. 
Let us here briefly discuss the decay products of the domain walls. After the collisions 
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of the domain walls, the energy stored in the walls is assumed to decay into the hidden 
massless U(l) gauge bosons. Then there is a relation between the energies of the hidden 
radiation and the gravitational waves: 

fl W ~ 2 (^^) ^, (26) 

where flhid is the density parameter of the hidden radiation. Therefore, the hidden radi- 
ation is always more abundant than the gravitational waves (see (|15p ). For the largest 
allowed value of A, we have flhid ~ 20fi gw . It might be possible to measure Q,ud indi- 
rectly by using the big bang nucleosynthesis (BBN), the cosmic microwave background, 
and large-scale structure data. Indeed, measurements of the light-element abundances 
give Vthidh 2 ^5 1.5 x 1CT 5 at 95% C.L. [36], which is translated to the bound on the gravi- 
tational waves, 

A 3 

n gw h 2 < 8 x 1(T 6 -=-. (27) 

Therefore, if the domain walls mainly decay into the hidden massless U(l) gauge bosons, 
the indirect measurement of Vt hid can set a bound on Q gw §. 

If SUSY is realized at a scale much higher than the weak scale, it may be difficult to 
directly produce SUSY particles in the accelerator experiments like LHC or ILC Then, it 
will be interesting if we can probe high-energy physics inaccessible to the accelerator ex- 
periments, by observing the gravitational waves. For m^^ ~ 10 4 — 10 6 GeV, the frequency 
of the gravitational waves is in the band of 10 2 Hz ^ / ~ 10 4 Hz, where the ground-based 
experiments are sensitive (see Fig. [I]). If the gravitino mass is heavier than 10 6 GeV, 
however, fo exceeds the frequencies to which the ongoing and planned experiments are 
sensitive. Even in this case, one may be able to indirectly observe the hidden radiation 
related to the gravitational waves by ff26l) . 

We have assumed that the constant in the superpotential, Wq, is the relevant source of 
the discrete R symmetry breaking. If the SU(N) gauge sector has couplings to the SUSY 
breaking sector, the gaugino mass rn\ generically receives contribution larger than the 
anomaly mediation. In that case the relation between fo and 1113/2 will be lost. However, 



8 It should be noted however that the presence of the hidden radiation is not an inevitable consequence. 
If the domain walls mainly decays into the visible particles, the bound (|27p is not applied. 
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the properties of the gravitational waves we have derived in the previous section are still 
valid if one replaces the gravitino mass by the (167r 2 /3iVg 2 )mA. 

Gaugino condensation is ubiquitous, if there are many gauge groups besides the SM 
ones in nature. The domain walls are generically formed when the gauginos condense. 
To avoid cosmological disaster, however, those domain walls must disappear. Since the 
discrete R symmetry is explicitly broken by the constant term wq in the superpotential, 
the domain walls are unstable and eventually decay into light particles. Interestingly, the 
gravitational waves are generated at the domain-wall collisions, and they contain valuable 
information on high-energy physics and the early universe. In this paper we have pointed 
out that the resultant gravitational waves can be a probe for measuring the gravitino 
mass, if wq is the relevant source of the discrete R symmetry breaking. 
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